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Abstract 

By parametrizing the action integral for the standard Schrodinger equation we present a derivation 
of the recently proposed method for quantizing a parametrized theory. The reformulation suggests a 
natural extension from conventional to nonlinear quantum mechanics. This generalization enables a 
unitary description of the quantum evolution for a broad class of constrained Hamiltonian systems 
with a nonlinear kinematic structure. In particular, the new theory is applicable to the quantization 
of cosmological models where a chosen gravitational degree of freedom acts as geometric time. This 
is demonstrated explicitly using three cosmological models: the Friedmann universe with a massless 
scalar field and Bianchi type I and IX models. Based on these investigations, the prospect of further 
developing the proposed quantization scheme in the context of quantum gravity is discussed. 

PACS numbers: 04.60.Ds, 04.20.Fy, 04.60.Kz, ll.10.Lm 



1 Introduction 

The Hamiltonian formalism and canonical quantization of Einstein's theory of gravitation was origi- 
nally developed by Dirac pQ, Arnowitt, Deser, Misner (ADM) [2], Wheeler DeWitt [I] and others. 
The Dirac- ADM action for general relativity is constructed by representing the spacetime metric us- 
ing the spatial metric components together with the lapse and shift functions that specify a time 
foliation. The latter set of functions enter into the action principle as Lagrangian multipliers ensur- 
ing the general covariance of the theory. Indeed, the 'already parametrized' form of the Dirac-ADM 
action requires that not all three-metric components and their conjugate momenta are to be freely 
specified on a spatial Cauchy surface. The Hamiltonian and momentum constraints must be satisfied 
initially and will be dynamically preserved A major challenge in quantizing general relativity is 
the treatment of these constraints. The straightforward application of the Dirac constraint quantiza- 
tion leads to a functional form of the Klein-Gordon equation, the Wheeler-DeWitt equation, and the 
associated difficulties in admitting a probabilistic interpretation. An alternative approach, the ADM 



1 



quantization method, attempts to circumvent these problems by isolating four out of the six spatial 
metric components (and their momenta) and regarding them as the evolution and coordinatization 
parameters, known as the embedding variables [51IH1II3, for the quantum evolution of the remaining 
'true' gravitational degrees of freedom. However different set of embedding variables can result in 
inequivalent quantizations. A further drawback of this approach is that the reduced Hamiltonian 
operator may involve a square root of an indefinite expression and therefore may not be Hermitian. 

The above pioneering work was followed by Kuchaf 's significant contributions to canonical quan- 
tum gravity in which Dirac quantization and the use of embedding variables are assimilated [7J |HJ |U] . 
Guided by the analogy with parametrized field theory, a canonical transformation is sought such that 
in the new set of variables the Hamiltonian and momentum constraints are linear in the conjugate 
momenta of the embedding variables |10| 111) . Remarkably, if this 'Kuchaf transformation' exists 
then Dirac and ADM quantization in fact coincide, with both methods leading to the same func- 
tional Schrodinger equation for the quantum evolution of the non-embedding variables. Furthermore, 
this approach holds out the possibility of constructing a framework for quantum geometrodynamics 
with the genuine Lie algebra of spacetime diffeomorphism. This was meticulously investigated by 
Isham and Kuchaf ^2]- A stringent criterion for Kuchaf s embedding variables is that they must be 
spacetime scalars, to guarantee the integrability of quantum evolution and its independence of time 
slicing Nevertheless, Kuchaf successfully demonstrated that in certain midisuperspace models, 
i.e. simple field-theoretical models with Killing symmetries, it is possible to find embedding variables 
explicitly that fulfill this requirement [HE]- Unfortunately, for some other important models, notably 
the 'scalar geons' |14H15lHHlll7j . finding the embedding variables has proven to be a formidable task. 

Recently, Hajicek and Kijowski developed a new approach in which the local existence of covariant 
gauge fixings is exploited to effect the decomposition of canonical variables on the constraint surface 
(defined by the Hamiltonian and momentum constraints) into embedding and dynamical variables 
|19U2Uj . The resulting 'Kuchaf decomposition' may form a basis for a complete Kuchaf transformation 
for a large class of geometries. The usefulness of this approach has been demonstrated by Hajicek and 
Kiefer in studying the quantum dynamics of a collapsing shell of null dust with spherical symmetry 
derived from the 'Louko- Whiting- Friedman model' [21] ■ However if a full theory of quantum gravity 
is to be constructed based on the covariant gauge fixing description then some obstructions will be 
encountered as have been noted in |2U| : A complete Kuchaf transformation becomes unattainable for 
geometries with Killing symmetries, which has also been observed by Torre [2^]. Besides, different 
gauges can lead to inequivalent quantum theories [21] and global gauge fixing may not exist at all 

[lElEE]. 

The purpose of this paper is to present an alternative approach to the quantization of a parametric 
theory, with a view to quantizing gravity, that will complement the above developments. While 
much insight has been gained from the quantization of a theory parametrized at the classical level, 
guidance for quantum gravity gleaned this way appears to be limited. Apart from technical difficulties 
and topological restrictions of the explicit constructions of Kuchaf 's embedding variables, there is a 
conceptual issue concerning these preferred variables. This is part of what has been attributed to 
the 'problem(s) of time' I26|. More specifically, the existence of a functional time relies on the 
canonical transformation for the geometro dynamical variables: (g a b,P ab ) {X A , Pa, Q r ,Pr) that splits 
the Hamiltonian and momentum constraints Ha(X b , Ps,q r ,p r ) into a 'kinematic part' 7i^(X B , Pb) 
and a 'dynamic part' Ti^(X B ,q r ,p r ), i.e. 

H A (X B ,P B ,q r , Pr ) = H A l (X B ,P B ) + H A >(X B ,q r ,p r ). (1) 

(Here a,b = 1,2,3; A,B = 0,1,2,3; r = 1,2.) A crucial assumption made by Kuchaf is that 
TL\(X b ,Pb) are linear in Pa so that the theory is 'deparametrizable' with X A and q r identified 
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as embedding and dynamical variables respectively. However an obstacle typical of midisuperspace 
models that prevents deparametrization is the nonlinear nature of Ti^(X B , Pb). Although in some 
cases, an ingenious choice of X A may well be the way out, this can be technically infeasible and might 
not be possible in principle. In these circumstances, we adopt the viewpoint that it may be possible 
to formulate a parametric quantum evolution of the dynamical variables q r and treat the remaining 
variables X A as unquantized kinematic variables without further restrictions. The new formulation 
sought should, nevertheless, reduce to Dirac quantization and is expected to be equivalent to certain 
other schemes up to a factor ordering |27| (in a semiclassical limit) if Kuchaf transformation exists 
that reduces H 1 §(X A , Pa) to simply Pa- 

As a first step towards the fully quantum field-theoretical formulation of our approach, a quantum 
mechanical description is adopted. In this simplified picture it is possible to focus on conceptual 
issues while establishing the essential methodology with a view to further generalization. Following 
a recent publication [2Sj where a new quantization scheme was proposed, we report in this paper 
how the proposed scheme can arise from a direct parametric description of quantum evolution as 
opposed to seeking to quantize a parametrized classical theory in a conventional manner. In section 
a parametric form of quantum evolution equations is introduced based on a parametrized action 
principle for the Schrodinger equation describing a non-relativistic particle. The advantage of this 
reformulation is that it opens up new avenues for generalized parametric quantum evolution beyond 
the description of a Newtonian particle. The generalized method is applied in section |3] to a simple 
Friedmann universe filled with a massless scalar field. Two further cosmological models, of Bianchi 
types I and IX, both involving gravity only are analyzed in sections fusing the nonlinear quantization 
scheme. Conclusions and discussions are made in section [3 where speculations on future work are 
discussed. Units in which c = h= 16irG = 1 are adopted throughout. 



2 Parametrized quantum evolution of a non-relativistic 
particle 

Consider a non-relativistic quantum particle of unit mass moving in an n-dimensional Riemannian 
manifold M with a time-dependent metric and coordinates q a , (1 < a < n) subject to a potential 
V(q a ,T). In terms of the Newtonian time T, the weight \ wavefunction ty(q a ,T) of the particle 
satisfies the Schrodinger equation 

i =h^f (2) 
dT K ' 

where the h = h(T) denotes the time-dependent Hamiltonian operator 

h(T) = -~A + V (3) 

with A denoting the Laplace-Beltrami operator on weig ht \ functions on M fTUIEK] . 
It has long been known that, in terms of the equal time inner product 

(*1,* 2 ) := U\^2 An q (4) 



of wavefunctions and the Schrodinger equation can be derived from the action integral 

(e.g. my- 

S Q [^,y*} = JdTM(y,(i&r-h)V} (5) 
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by freely varying ^ and its complex conjugate ^f*. 1 However the potential significance of this formula- 
tion in the context of constrained Hamiltonian systems does not appear to have been fully recognized 
in the literature and will be explored in this paper. In what follows an action integral equivalent to 
(|5|) will be constructed that generates quantum evolution equations with respect to a general time 
coordinate, by parametrizing the Newtonian time and turning it into a constrained classical variable. 

We anticipate the parametrized T and its momentum to couple with a reformulated Schrodinger 
equation in a manner similar to the interaction between classical and quantum variables. Of course, 
the 'classical variables' arising from this procedure do not carry additional physical degrees of freedom 
due to the constraining relations accompanying the parametrized theory. Therefore these 'constrained 
classical variables' are really kinematic (embedding) variables, rather than genuine dynamical vari- 
ables, forging a link between the physical and geometrical descriptions of the systems |28j . Nonetheless 
the expected mathematical similarity with the semiclassical theory suggests that the action integral of 
the parametrized quantum theory is a sum of its kinematic and dynamic parts. While the kinematic 
part depends on 'classical variables' only the dynamic part will involve both classical and quan- 
tum variables derived from the original action ©, whose solutions have the scaling invariance under 
^ — > for any complex constant A, even though |A| 7^ 1 [3131 • However this scaling invariance will 
apparently be violated due to the presence of the kinematic part of the action that is independent of 
the quantum variables. A way to avoid this problem is to eliminate the scaling invariance by using, 
instead of (JSJ) , the modified action integral [3UJ 1^ : 

S Q [^,^*,a] = JdT{$t{ty,id T ^) - (V,hifr) + a ((*,*) - 1)} (6) 

with a = a(T). With respect to this Lagrangian multiplier, the variation of © enforces the normal- 
ization condition 

(*,*) = 1 (7) 

for the wave function whereas the variation with respect to \E' and \&* yields the Schrodinger 
equation of the form 

i^ = (h-a)9. (8) 

Despite that a can be specified arbitrarily, different choices of it correspond to wavefunctions unitarily 
related by a phase depending only on time and therefore considered as physically equivalent. 

The action ((HJ) will now be brought to a parametric form. This procedure starts by expressing T 
as an arbitrarily chosen function of a parameter time t to yield 

5 Q [*,^*,a]= |di{^,0 t *)-r[(^,^*)-a((*,*)-l)]} (9) 

where the over dot denotes dt- Note that in this expression = *$>(q a ,t) := fy(q a ,T(t)), a = a(t) := 
a(T(t)) and h = h(t) := h(T(t)) (with slight abuse of notation). At this stage T{t) is merely an 
externally supplied function. It can, however, be included as part of an enlarged set of dependent 
variables augmented with a new pair of Lagrangian multiplier and constraint. To this end introduce 
the variable n(t) and use it to define the quantity 

Hq :=n+(*, h *)- a ((*,*)- 1) (10) 



-'The subscript 'Q' indicates 'quantum'. 
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so that, if TCq vanishes, II equals the coefficient of T in the action thereby serving as the momentum of 
T(t). By replacing the coefficient of T with II in the Lagrangian for Q and adjoining this Lagrangian 
with the term —NTCq using a Lagrangian multiplier N = N(t), we are led to the following action 

S Q [y,V*,a,T,IL,N} = Jdt{R{y,id t *) + nf-NH Q }. (11) 

This action takes a parametric form and is equivalent to © upon variation with respect to all 
components in the extended set of variables {vf, VP*, a, T, IT, N}. It is not hard to further express this 
action in a complete canonical form by defining the conjugate momenta of \P and VP*. This form is 
however not required in our present analysis. 

It is instructive to split TCq into its 'kinematic' part 

Hi = 11 (12) 

and 'dynamic' part 

W§ = (*,fe*)-a((* J *)-l) (13) 

so that 

h q = h% + h%. (14) 

Accordingly, the structure of parametrized action (|11|) can be made more transparent with a similar 
split: 

Sq = S$ + S q > (15) 

in terms of the kinematic part 

S%[T,I1,N} = Jdt{nt-NH%} (16) 

and the dynamic part 

S%[%**,a,T,Il,N] = Jdt{M(V,id t y)-NH%}. (17) 

The parametrized quantum evolution equations can be explicitly derived from the action (|11|) under 
variations with respect to and its conjugate, a, T, H and N to be 

idt® = N(h - a) * (18) 



(*,*) = 1 (19) 
t = (20) 

n = -n 9 ^ (2i) 

n Q = (22) 

respectively. The structure of these equations indeed resembles that of a semiclassical theory. It is 
evident that the kinematic action Sq involves only 'classical' variables and takes a canonical form. 
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This makes it possible to perform a canonical transformation, say, from (T, II) to (T' , IT) by preserving 
the form of <Sq [T, IT, N] = Sq[T', IT, N]. It follows immediately that with these transformed variables, 
the action Sq^, \E'*, a, T 1 , II', N] will yield equivalent quantum evolution equations of the same form 
as H18p - (j2'2|) with the substitutions T — > T' and II — ^ II'. This suggests the set of parametric quantum 
evolution equations is somehow more flexible in choosing a 'time' variable than that of the Schrodinger 
equation @ that assumes the existence of a preferred choice of time. A further advantage of the 
present approach is the implication that the action (jllj) for a wider class of the kinematic actions 
with a general = Hq(T, II) in place of (HJ) may play a fundamental role in the quantization of a 
parametric theory where no preferred (Newtonian) time can be identified. In this case the parametric 
quantum evolution system (|18j) - ()22|) becomes nonlinear and cannot be reduced to the Schrodinger 
equation ((2j). The following sections will provide illustrative examples for this generalization based 
on minisuperspace models. 



3 Friedmann universe with a quantized massless scalar 
field 

The theoretical framework for the quantization scheme developed in the previous section will now 
be employed in the quantization of a simple Friedmann universe filled with a massless scalar field. 
The nonlinear quantization of a similar model with a massive scalar field was addressed in |2H] in 
terms of a nonlinear integro-partial differential system derived from (|18f) — (|22|) (for a = 0), which was 
reduced to a system of nonlinear ordinary differential equations of infinite dimensions. The latter was 
further truncated for numerical simulation. Not surprisingly, the reduced system can be shown to 
be non-dissipative since there exists an underlying variational principle as developed in the previous 
section. 

The present Friedmann model with a massless scalar field can be investigated analytically. The 
classical Lagrangian for this model may be found from |32| and |2H] by dropping the mass term to 
obtain 

L((j),4>,R,R,N) = -—R 2 + —4> 2 + 6NKR. (23) 

Here (ft = <j)(t) is the scalar field <j> which has been conveniently re-scaled. The positive variables R = 
R(t) and ./V = N(t) are the scale factor and lapse function respectively appearing in the Robertson- 
Walker metric 

ds 2 = -N 2 dt 2 + R 2 d£ 2 (24) 

with d£ 2 denoting the squared line element on the homogeneous and isotropic 3-space, which may be 
closed (K = 1), flat (K = 0) or 'open' (K = — 1). The momenta conjugate to R and </> are 

^ dL 12R- 

n := — r = t^R 25 

dR N v ; 

p: =5r— *■ (26) 

It follows that the corresponding Hamiltonian H((fi,p,R,H,N) := pep + UR — L takes the form 
H = NH where 

n<t,P,R,m = ~ + ^-eKR. (27) 



6 



The canonical equations of classical motion can be derived with arbitrary N(t). The gauge condition 
N = 12i? 3 is chosen so that these equations take a simpler form as follows 



= P (28) 
p = (29) 
R = -R 2 U (30) 

(31) 



n = -^! + ^ + 72^ 

2 2R 



subject to the (Hamiltonian) constraint 



n 2 v 2 

+ ^-^KR = 0. (32) 



24R 24R 3 

Clearly, from Q29|) the momentum p is a constant of motion, with which (|28|) is readily integrated. 
With a suitable choice of the origin of the coordinate time t, one obtains 

R 2 = gl (33) 
6(1 + KeW) V 7 

for — oo < t < oo with X = 0, 1 and — oo < t < with = —1. In terms of this expression (|3fl|) can 
be easily solved for II. 

The quantization of this model requires a degree of freedom to be quantized. We choose this to 
be the scalar field (j), thereby regarding the scale factor as the 'geometric time'. As per discussions in 
section |2] we split TL by 

n = n K + n D (34) 



where 



To proceed the operator 



is constructed by substituting p — > p := — i^ into W D . This operator is Hermitian with respect to 
the inner product 

poo 

<#i,tf 2 > := / *i*2d0 (38) 



of any two wavefunctions ^i(</>, t) and ^(^j^)- 

In the gauge N = 12i? 3 , the following nonlinear evolution equations for the normalized wavefunc- 
tion ^(4>,t) and kinematic variables R(t),H(t) then arise from (fT5 |) -(|2*2* |) : 



1 d 2 ^ 



i- 



dt " 2 d6 2 



(39) 



R = -R 2 U 



(40) 



n = - 



RU 2 3 (p 2 ) 



+ 72KR 3 



(41) 



2 2R 



subject to the constraint 



n 2 (p 2 ) 



- 6KR = 



(42) 



24R 24R 3 



where 



■oo 



(43) 



and the phase condition a = —6KR has been used for simplicity. These equations are similar in 
structure to those obtained by Kim in analyzing the Friedmann universe filled with a massless scalar 
field in the context of semiclassical gravity |33j . 

The wave equation ()39|) in fact decouples from the others and takes the same form as the 1- 
dimensional Schrodinger equation for a free non-relativistic particle. Hence (]3 2 ) is preserved under 
evolution. It follows that the solution for R can be generated from the classical solution ()33|) by 
replacing \p\ with \J (p 2 ). Like the classical case, there is no singularity avoidance and R can become 
arbitrarily close to zero for all possibilities of K. Unlike the classical case, however, it is possible 
to envisage a wave packet with zero mean scalar field and momentum ((</>) = (p) = 0) but nonzero 
deviations {{4> 2 ) > and (p 2 ) > 0). In this case, the evolution of the Friedmann universe can be 
thought of as being 'purely quantum-driven'. 



4 Nonlinear quantum cosmology of the Bianchi type 



The foregoing Friedmann universe with a scalar field was quantized while maintaining the classical 
nature of the scale factor which carries the gravitational degree of freedom of the universe. In this 
section, two vacuum cosmological models, of Bianchi types I and IX, will be examined. Both belong 
to a wider class of homogeneous but anisotropic spacetimes |34j . In Misner's seminal paper [221, the 
classical dynamics and quantization of Bianchi I and IX universes in the context of ADM quantization 
has been studied in detail. Below we will give derivations of the quantum evolution equations for 
these two cosmological models and comment on the consequences. The emphasis will be placed on 
the nonlinearity inherent in the derived quantum systems. 

Just as in the Friedmann universe, the geometrical description of the Bianchi I and IX models 
involves the lapse function N(t) and scale factor R(t). Two additional functions f3+(t) and @-(t) 
are used, to allow for the dynamical anisotropy of the spacial hypersurface at any coordinate time t. 
Following Misner's notation we consider a spacetime metric of the form 



for i,j = 1,2,3, where cr, denote some basis 1-forms of the spatial hypersurface and (3 is a traceless 
matrix with elements given by 



The constant time spacial hypersurface for the Bianchi I universe is open and has an Abelian homo- 
geneity group. Thus Oi are chosen so that dai = and can be simply expressed as a\ = dx, 02 = 




(44) 




(45) 
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dy, C3 = dz, for — oo < x,y, z < oo. The homogeneity group of the closed spacial hypersurface for the 
Bianchi IX universe is that of S 3 . Therefore the basis 1-forms Oi are chosen to satisfy the structure 
equations 



1 

2 

and can be expressed explicitly as 



doi = - e ijk Oj A <r k (46) 



o"i = sin ipdi9 — cos ip sin fidcp, 02 = cos ipd'd + sin ip sin fidip, 03 = — dip — cos (47) 

for < ■& < ir, < ip < 2ir and < ip < 4tt. 

The classical dynamics of the models is generated by the Dirac-ADM action subject to the homo- 
geneity condition for all metric functions. Up to an overall factor, the resulting Lagrangian takes the 
form 

L(J3±,$±,R,R,N) = -^R 2 + ^ ($1 + (3 2 ) + NR 3 R (48) 



N N 

where 1Z denotes the scalar curvature on the spatial hypersurface. It can be conveniently expressed 
as 

n =^- V ) ( 49 ) 

where 

V = 1 (50) 

for the Bianchi I model (such that 1Z = 0) and 
V = V(p+,/3-) = 

1 + I e -8/3+ + I 4(/3 ++v / 3/3_) , I 4(/3+- v / 3/3_) _ 2 4/3+ _ 2 2 (f3++V3(3-) _ I -2(J3+-y/3/3-) (kiN 

3 3 3 3 3 3 v ' 
for the Bianchi IX model. By introducing the conjugate momenta 

n := *l = (52) 

OR N { J 

p±:= m = ^r p± (53) 

we can construct the Hamiltonian H(/3±,p±, R,IL, N) := HR + p + j3 + +p_/3_ — L = NTi where 

n(P ± , P± ,R,U) = - J^L + ^ (p 2 + +p 2 ) + ^ (V- 1). (54) 

From (|54|) the canonical equations of motion can be derived and take a simpler form using the choice 
N = 12R 3 as in the Friedmann case. In this gauge, denoting d± := dg,, one obtains: 

P± = P± (55) 

p± = -im i d±V (56) 

R = -ILR 2 (57) 
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n = ^ + ^ {p% + pl) - 18R 3 (V - 1) (58) 
and the Hamiltonian constraint 

_In 2 i? 2 + \{p\ + pi) + i8i? 4 (y - 1) = o. (59) 

The analogy between the Bianchi and Friedmann models should be apparent. Specifically by com- 
paring the Lagrangian in ()48j) with that in (|23|) we see that the kinetic terms in both cases share the 
same structure with the functions (3± resembling the scalar field <p. 

For the Bianchi I model with V = 1, both p + and p_ are constants of motion. By introducing 



P-=^P 2 ++P- (60) 

the solution for R may be generated from the scale factor in the Friedmann model with zero curvature, 
i.e. K = 0, given in (|33|) . This yields 

R 2 = V -e 2 ^ (61) 

for — oo < t < oo, from which II can be evaluated using (|57[). The scale factor (|61|) together with 
N = 12R? and arbitrary constants p± is equivalent to the Kasner solution for the Bianchi I universe 

m 

The quantization of the Bianchi models can be performed once the dynamical variables to be 
quantized are decided. Guided by the analogous roles played by (3± and 4>i we choose to turn f}± 
into quantum variables and introduce their wavefunctions \&(/3±, t). The inner product of two such 
wavefunctions ^ can be defined to be 

rco /'OO 

<tfi,tf 2 > := / / *I* 2 d/? + d/3_. (62) 



— ooJ — oo 



We are concerned only with normalized wavefunctions ^ satisfying (^Sf, \P) = 1 so that (6) = (\E , ,o\E f ) 
for any operator o. 

As before, the scale factor R will act as the classical time that evolves with respect to the coordinate 
time t. The coupling between these quantum and classical variables follow from ()18j) ~ (j22|) applied in 
the current setting. In accordance with procedures leading to ()39j) - (|42|) . we perform the 'kinematic- 
dynamic' split 

n = n K + n D (63) 



where 



""-ss (64) 



and construct the operator 



nD =2^^ +p2 ^ +3 4 {y - i) (65) 

1 



h = -^(dl + di)^ + 18R 4 (V -I) (66) 

from HP with the substitutions p± — > p± := —id±. 

By using the conditions N = 12i? 3 and a = we arrive at the following nonlinear evolution 
equations 

id t qj = --(d% + d 2 ) ^ + 18R 4 (V -1)^ (67) 
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R = -TIE 2 



(68) 



n = 



(69) 



subject to the constraint 



|n 2 i? 2 + 1 (f + + f_ > + i 8j R 4 ((y) - 1) = o 



(70) 



for *(/?±,t), i2(t) and H(t). 

For the quantized Bianchi I model (V = 1), the decoupled wave equation (|6*7j) has the same 
structure as the Schrodinger equation for a free non-relativistic particle in 2-dimensions as expected. 
Hence {p\ + p 2 _) is preserved under evolution and the solution for R can be generated from the 



classical solution (|6*Tj) by replacing p with J {p 2 + +p 2 _). Clearly as t — > — oo the scale factor tends 
to zero, i.e. R — > 0, where the spacetime becomes singular. Using the above properties, one may 
demonstrate explicitly how this universe may be driven by quantized geometry described by a wave 
packet with zero means of the dynamical gravitational fields and their momenta ((/3±) = (p±) = 0) 
but with nonzero deviations ((/3+ +/??.) > and (p\ +pL) > 0). 

For the quantized Bianchi IX model where V is given in (|51|) . equations l|67 |) -(|7U| l constitute a 
system of fully coupled nonlinear integro-partial differential equations. In the absence of explicit 
solutions it is not possible to analyze the dynamical behaviour of this model in detail. Nonetheless, 
it is worth noting two generic properties that one may conclude from this system. First, it can be 
shown that V is nonnegative with its minimum V = taking place at f5± = 0. (In fact V ~ 8f3 2 _ + 8(3 2 _ 
for \/3±\ <C 1. For further details of V, see It follows that the term V — 1 in TC can become 

negative which could be problematic for ADM quantization due to a square root taking procedure. 
The proposed quantization scheme is free from this problem, as exemplified in ()67|) - ()7U|) where no 
square roots are involved. Secondly, although the scale factor is treated as a classical variable it 
is not used as an evolution parameter. As such, the quantum evolution will not come to a halt at 
the 'maximum hypersurface' where R = classically. By comparison, the pathology associated with 
time integration across such a hypersurface could be suffered by both ADM quantization and Dirac 
quantization using Kuchaf's embedding variables 



By virtue of its action principle, we have reformulated the Schrodinger equation in non-relativistic 
quantum mechanics into a parametric form. In this formalism the Newtonian time acts as a kinematic 
variable that evolves with respect to a general time coordinate while coupling to a quantum wavefunc- 
tion in a semiclassical fashion. A gratifying feature of this description is that the parametrized action 
decomposes into a kinematic part and dynamic part. The former part of the action depends only on 
classical variables and takes a canonical form. Therefore, by preserving this form, one may perform a 
canonical transformation that redefines the kinematic time and its conjugate momentum. This sug- 
gests a natural extension of the proposed quantization method by accommodating a general expression 
of the kinematic time and its conjugate momentum that enter into the kinematic action. A quan- 
tum description obtained from this generalization is no longer reducible to the standard Schrodinger 
equation. Instead, it leads to a theory for nonlinear quantum evolution previously proposed based on 




5 Conclusions and discussions 



n 



heuristic arguments |28j . In spite of the breakdown of the principle of superposition in time evolution, 
the resulting theory belongs to a class of nonlinear quantum theories admitting both probabilistic 
interpretation and a Hilbert space with Euclidean norm |37l 138) . 

These properties have indeed been demonstrated in the quantization of three types of cosmological 
models in this paper. The first model we considered was the Friedmann universe with a massless scalar 
field. The model consists of both gravitational and matter degrees of freedom, carried by the scale 
factor and scalar field respectively. Although classically these two variables take the same status as 
dynamical variables, given the physical origin of this model, we believe that it is matter that should be 
quantized thereby making the scale factor the intrinsic geometric time. In doing so we have not treated 
the scale factor merely as an evolution parameter as in certain alternative quantization methods. One 
obvious benefit of allowing the intrinsic geometric time to evolve as a classical variable is so that the 
theoretical description is not confined to scenarios with monotonically ascending geometric time. This 
avoids the diverse choice of time variables for different curvature parameters if the same cosmological 
model was quantized using the ADM method |32| . It might be perceived at this point that our 
approach only applies to the quantization of matter in classical spacetime. To explicitly demonstrate 
it is not the case, we then considered the Bianchi I and IX cosmological models describing the evolution 
of an 'empty' universe. No matter fields are present there and one must decide which gravitational 
degrees of freedom should be quantized. In section |1] the two metric functions representing the 
anisotropy of the universe have been chosen for this purpose. Arguably this is the most appropriate 
choice due to the physical interpretation of these two functions as the gravitational 'waves' in the 
cosmological model. These considerations suggest that in extending the proposed quantization to 
quantum gravity, a primary concern is to determine which two gravitational field components are 
to be quantized so that the remaining four fields will act as classical kinematic variables. Unlike 
Kuchaf's embedding variables, however, our approach does not require that the momenta of the 
kinematic variables to be resolved in the Hamiltonian and momentum constraints. In contrast, it 
may be speculated that if gravity were to be quantized based on the proposed quantization then the 
corresponding action integral should contain a kinematic part of the form 

S K [X A ,P A ,N A ] = JdtJ d 3 x{p A X A - N A H\{X B ,P B )} (71) 

with A,B = 0, 1, 2, 3 where N A are Lagrangian multipliers proportional to the lapse and shift func- 
tions, X A and Pa denote the gravitational fields and their momenta to act as classical kinematic 
variables and TL\{X b : Pb) are some expressions, not necessarily linear in Pa, to be determined. If 
such a split can be found uniquely in the full theory, then the classification of the geometrodynamical 
variables into kinematic and dynamic sets will follow naturally. In contrast, the classifications of vari- 
ables in the truncated cosmological models in sections El and 0] have been partly guided by physical 
intuition. 

A main thrust for this work has been the preservation of unitarity in quantizing physical systems 
with generic nonlinear kinematical structure. Of course, there exist other model-dependent strategies 
for unitary evolution within the linear quantum formulation, e.g. |21j . In a broader context, even 
the principle of unitarity can be debated by relaxing the probabilistic interpretation. Some have 
suggested that it is at best an approximate concept [331 EH- However, a more conclusive argument 
on the probabilistic interpretation of quantum gravity does not appear to be available at present. 
The non-unitary Wheeler-DeWitt equation is thought to be capable of describing a tunnelling |1J 
or topology-changing [121 wavefunction in certain models of the early universe. To what extent these 
theoretical scenarios can be accommodated or precluded by our approach would form a subject for 
future investigation. 
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Although the proposed framework is somewhat reminiscent of semiclassical gravity original devel- 
oped by M0ller and Rosenfeld [35] as an exact theory, where all gravitational degrees of freedom are 
classical, it is worth stressing that we seek to quantize two gravitational field components. On the 
other hand, the remaining kinematic variables in our approach are classical by construction and there- 
fore do not arise from a decoherence process required by Kiefer's approach to semiclassical gravity as 
an approximation to quantum gravity |4(J| . 

Before further potential implications of this work on quantum gravity can be contemplated se- 
riously, several pressing issues ought to be addressed first, regarding the extension of the current 
quantum mechanical approach to a field-theoretical description. In this respect, it would be prof- 
itable to investigate the quantum dynamics of systems such as scalar fields in curved spacetime and 
gravitational wave spacetimes [IHj • Progress made along these lines will pave the way for the analysis 
of a number of midisuperspace models where Dirac quantization using Kuchaf 's method has not been 
successfully applied. Results from research described above will be reported elsewhere. 
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